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ABSTRACT 



^ ' We present a parsec-scale image of the OH maser in the nucleus of the active 

■ galaxy IIIZw35, made using the Very Long Baseline Array at a wavelength of 

18 cm. We detected two distinct components, with a projected separation of 
50 pc (for D = 110 Mpc) and a separation in Doppler velocity of 70kms~^, which 
contain 50% of the total maser flux. Velocity gradients within these components 
could indicate rotation of clouds with binding mass densities of ~ 7000 Mq pc~^ 
and total masses > 5 x 10^ Mq. Emission in the 1665-MHz OH line is roughly 
coincident in position with that in the 1667-MHz line, although the lines peak 
at different Doppler velocities. We detected no A = 18 cm continuum emission; 
our upper limit implies a peak apparent optical depth ItohI > 3.4, assuming the 
maser is an unsaturated amplifier of continuum radiation. 



Subject headings: masers — galaxies: active — galaxies: nuclei — galaxies: 
individual (IIIZw35) 
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1. Introduction 

IIIZw35 is a pair of galaxies catalogued in Zwicky (1971) as a 16th-magnitude compact 
blue double. Its systemic velocity, based on observations of CO emission ( |lVlirabel~SZ 



Sanders 19871 ), is 8289 km s ^ (heliocentric, optical definition), corresponding to a distance 



of 110 Mpc {Hq = 75kms ^ Mpc ^; Imas is equivalent to a projected separation of 0.53 pc). 
Optical and IR images ([Chapman et al. 1990|) show two components separated by 9" along 



a position angle of 20°. In optical images, the northern galaxy is brighter and redder than 
the southern, and emits the majority of its energy in the far infrared; IRAS measured its 
10 — 120 fim luminosity to be 3.5 x 10^^ Lq, with a spectral energy distribution consistent 
with thermal emission from dust at 65 K. Chapman et al. (1990) classify the northern 
galaxy as either a LINER or Seyfert 2, based on optical emission line ratios. 

The nucleus of the northern component of IIIZw35 is the site of strong OH maser 
emission. The total luminosity in the 1667-MHz transition, assuming isotropic emission, 
is 530 Lq, with a ratio of flux density in the 1667 and 1665-MHz spectral peaks of 9:1 
( ^taveley-Smith et al. 1987|) . The 1667-MHz spectrum has two peaks at approximately 
±35kms~^ from the systemic velocity. The maser extends in Doppler velocity over about 
lOOkms"^, but there is also faint (< lOmJy) blueshifted emission extending more than 
500kms~^ from Vgys- Baan, Haschick, & Henkel (1989) interpret the blueshifted wings in 
this and other OH maser galaxies to be due to molecular outflows. 

Chapman et al. (1990) imaged the A = 18 cm continuum emission from IIIZw35 with 
the VLA and MERLIN and measured a flux density of 30 ± lOmJy and an angular size of 
< 300 mas (Tg > 10'^ K). The centroid of the maser emission was observed to coincide with 
that of the 18 cm continuum to within 0.3". Assuming the maser acts as an unsaturated 
amplifier and completely covers the continuum, they inferred an apparent peak OH optical 
depth ItohI ~ Montgomery & Cohen (1992) mapped the OH maser emission with 
MERLIN and observed an extended (~ 80 mas) structure with an apparent velocity gradient 
of 2.5kms~^pc~^ along a position angle of —20°. They interpreted the gradient to be the 
signature of a rotating molecular disk in the nucleus of IIIZw35. 
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2. Observations 

We observed IIIZw35 on 25 March 1995 at A = 18 cm with the NRAO Very Long 
Baseline Array (VLBA)|^ in spectral-line mode. Only the stations at Fort Davis, Kitt 
Peak, Los Alamos, North Liberty, Owens Valley and Pie Town provided useful data. The 
baseline fringe spacings ranged from 17 to 200 mas. The data were recorded in both 
left and right circular polarization, in two 8-MHz passbands centered on f = 8438 and 
7000 km s~^. The former covers both the 1667-MHz and the 1665-MHz transitions of 
OH, while the latter was tuned to lie outside the known range of maser emission. The 
data were correlated at NRAO; each passband was divided into 512 spectral channels, 
yielding a Doppler velocity spacing of 2.8kms~^. All further processing was undertaken 
with AIPS at the CfA. After initial calibration of the amplitudes, fringe rates, and 
delays, the visibilities were phase-referenced to the channel of the peak line emission, at 
8245 km s~^. Analysis of residual fringe rates established the position of the maser peak: 
a^ggo = 01^41'^47?925 ± 0?007, 61950 = 16°51'05'.'8 ± 0'.'2. An error in this reference position 
during processing resulted in a systematic error in the positions of spectral features that 
scales linearly with frequency offset from the maser peak. The magnitude of this error is 
1.7masMHz~^; the resulting systematic errors are < 1 mas in the 1667-MHz line profile, and 
< 4 mas in the 1665-MHz line profile. The 1-cr noise level in a uniformly-weighted spectral 
channel image was 3.8mJy/beam. Spectral channels in the lower-frequency passband were 
averaged to produce continuum visibilities. The l-o" noise level in a uniformly-weighted 
continuum image was 0.18 mJy/beam. 

3. Results 

Figure 1 shows the distribution of 1667-MHz OH maser emission in the nucleus of 
IIIZw35. We detected emission in two clumps, with a projected separation of 90 mas (50 pc) 
along an axis with a position angle of —10°. The spectral peak of the maser emission is in 
the southern clump. Figure 2 shows expanded views of the northern and southern clumps. 
The positions of the 1665-MHz line peaks are also plotted. In the single-channel images, 
the maser features were unresolved; we place an upper limit of 17 mas on their angular 
sizes, which implies a peak brightness temperature Tqh > 3 x 10® K. Figure 3 compares the 
imaged flux density in both clumps with the single-dish spectrum of Staveley-Smith et al. 
(1987). Note that ~ 50% of the total maser flux density was resolved in our VLB A images. 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement with Associated Universities, Inc. 
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Sensitivity limitations did not permit us to image the full range of bluesliifted emission, 
which extends ~ 500kms~^ from the systemic velocity at a level of less than lOmJy ( |Baan, 
Haschick, fc Henkel 1989| ). 

The northern clump has a projected extent of ~ 10 pc (Figure 2a). It consists of 
multiple spectral features, over a velocity range of llOkms"^ (Figure 3). We observed 
an east-to-west velocity gradient of llkms~^pc~^ over ~ 5pc in the velocity range 
8290 <v < 8340 km s"^ (Figure 1, mset). Between 8230 and 8290kms"\ we did not detect 
a gradient; all of the emission in this range lies in a compact (< 2pc) cluster at the SE 
end of the clump (Figure 2a). A recent high-resolution intercontinental VLBI observation 
revealed a north-south gradient in this cluster of 30kms^^pc^^ (P. Diamond, private 
communication). Lower sensitivity and angular resolution prevented us from detecting a 
gradient of this magnitude. 

The brightest emission is in the southern clump. The imaged line profile is well fit by a 
Gaussian function with a peak of 130 mJy at 8245 km s~^ and a FWHM of 30kms~^, with 
line wings extending ±50kms~^ from the peak at a level of < 50mJy. Most of the emission 
in the southern clump is confined to a region approximately 3pc across (Figure 2b); 
however, emission in the line wings extends ~ 5 pc SE and NW of the maser peak. Evidence 
for a velocity gradient is ambiguous within the compact part of the southern clump, but 
emission from the line wings suggests a velocity gradient of about lOkms"^ pc~^, increasing 
along a position angle of 110°. This gradient is of roughly the same magnitude as that in 
the northern clump, but in the opposite direction (Figure 1). 

We imaged 1665-MHz maser emission in both the northern and southern clumps. The 
southern clump has two spectral peaks, at 8232 and 8258 km s~^, both with fiux densities 
of ~ lOmJy. These two peaks straddle the 1667-MHz peak in Doppler velocity (Figure 3). 
Both of the southern 1665-MHz features are coincident with the compact component of 
the 1667-MHz clump within the position errors (Figure 2b). In the northern clump, we 
detected one 1665-MHz feature, also with a peak fiux density of ~ 10 mJy. The velocity of 
this feature is 8353 km s~^, which is near the red edge of the 1667-MHz profile. Moreover, 
while the velocity of the 1667-MHz emission increases to the west in the northern clump, 
this 1665-MHz feature lies at its eastern end (Figure 2a). It appears that the two OH maser 
clumps differ in their relative distributions of 1667 and 1665-MHz emission, both in angle 
and in velocity. Note that the features we attribute to the 1665-MHz transition could be 
redshifted (~ 8600 km s~^) 1667-MHz features. Single-dish spectra {e.g., [Mirabel fc SandeTs 
|1987D do show faint emission filling the gap between the main 1667-MHz profile and the 
putative 1665-MHz profile. However, the near angular and Doppler velocity coincidence 
of emission in the two profiles suggests we are detecting the 1665-MHz transition in both 
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clumps, rather than highly redshifted 1667-MHz features. 

We did not detect any 18 cm continuum emission in the nucleus of IIIZw35 with 
the VLBA. We imaged various subsets of the data with a variety of weighting schemes 
to confirm this non-detection. We place a 5a upper limit on the continuum flux density 
of ImJy/beam. The visibility amplitudes, averaged over the entire observation, were 
consistent with a = ImJy noise on all baselines (180-2200 km). 



4. Discussion 

The concentration of maser emission into two compact clumps is likely a result of 
inhomogeneity in the distribution and excitation of OH in the inner 100 pc of IIIZw35, and 
not to compact structure in the background continuum. Chapman et al. (1990) observed 
a A = 18 cm continuum source with a flux density of ~ 30 mJy with MERLIN. The 
non-detection of the continuum on our shortest VLBA baselines places a lower limit on the 
angular size of a 30mJy Gaussian source of 170 mas, or an upper limit on the continuum 
brightness temperature, Tc < 10^ K. This suggests that the 18 cm nuclear continuum 
structure is smooth on the angular scale of the maser clump separation. It is also possible 
that filamentary compact continuum sources of < 1 mJy are being amplified by the northern 
and southern maser clumps. However, the components of the northern clump show a 
clear gradient of velocity with position, suggesting that they sample a coherent, possibly 
rotating, region of molecular gas (Figure 1, inset). 50% of the total maser flux density was 
resolved in our VLBA images (Figure 3), which may indicate an extended halo component 
of maser emission. Smoothly distributed molecular gas could act as a low-gain amplifier 
of an extended continuum source. Alternatively, there may exist a population of compact 
maser clumps, each with a flux density below our detection limit. 

If the velocity gradient of llkms~^pc^^ that we observe over 5pc in the northern 
clump (Figure 1; inset) is due to rotation of a homogeneous, gravitationally bound cloud, 
the implied mass density is ~ 7000 M0 pc-3 (total enclosed mass ^ 5 X 10^ M©). The 
gradient suggested in the line wings of the southern clump could be produced by rotation 
of a cloud of similar density. A difficulty in interpreting gradients in spectral images arises 
when two or more distinct components are close in angle (less than a beamwidth) and 
velocity (less than a linewidth). Blending of the individual components can simulate a 



smooth velocity gradient in an extended feature. The MERLIN observations of [Montgomery 



&; Cohen 1992| had a beamwidth of 280 mas and a channel spacing of 14kms ^; at this 



resolution, spectral blending between the two distinct maser clumps appeared as a gradient. 
True gradients can only be indentified with adequate resolution. With this caution in mind. 
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we note that the S-shaped structure in the northern maser clump (Figure 2a) could be 
produced by four or more blended, unresolved features. Even so, these features would still 
display a correlation between velocity and position. 

Baan, Haschick, & Henkel (1989) noted that the apparent luminosity of extragalactic 
OH masers scales approximately as L"^-^^. They interpret this correlation to imply that 
pumping by nuclear FIR radiation produces a column density of inverted OH that is 
proportional to LpiR, which then acts as a low-gain amplifier of the 18 cm radio continuum, 
itself proportional to LpiR ( [Norris, Allen, fc Roche 1988|) . The apparent optical depth 
estimated for the OH maser in IIIZw35, |toh| ~ 1-7, was based on the assumption that the 



maser completely covers the background continuum ( [Chapman et al. 1990] ). Our observation 



indicates that the true covering factor is significantly less than unity. We place an upper 
limit on the brightness temperature of the 18 cm background continuum, Tc ^ 10^ K, and a 
lower limit on the peak brightness temperature of the maser, Tqh > 3 x 10^ K. Assuming 
the maser is unsaturated, we derive a lower limit on the apparent optical depth, ItohI > 3.4. 

If we assume that the maser is an unsaturated amplifier of background continuum 
emission, that the excitation temperatures of the 1667 and 1665-MHz transitions are equal, 
and an LTE optical depth ratio ti667/ti665 = 1-8, then the ratio of 1667 to 1665-MHz line 
peaks in single-dish spectra of 9:1 implies |ti667| ~ 5 ([Henkel fc Wilson 1990|) . However, the 



two transitions do not peak at exactly the same velocities, and the resulting variation in line 
ratios suggests 3 < \tiqqj\ < 5 across the emission profile. Randell et al. (1995) compared 
the observed main-line and satellite-line ratios with models of extragalactic OH emission 
regions, and found that the OH line ratios are sensitive to a number of properties of the 
molecular gas, including density, dust temperature and the presence of velocity gradients. 
The observed range of 1667 to 1665-MHz line ratios may be due to variations of these 
quantities within the maser clumps. Our lower limit on the peak apparent optical depth, 
kieerl > 3.4, is consistent with the range of optical depths implied by the mainline flux 
ratios. It therefore appears that low-gain models of extragalactic OH masers {e.g. [Baan 
Haschick, fc Henkel 1989| ) do not apply in IIIZw35. 



One possible explanation for the angular structure of this OH maser is that the clumps 
trace edge-on segments of molecular shocks. Shocks are expected to be associated with 
energetic molecular outflows, such as those inferred from the blueshifted line wings in 
several OH maser galaxies, including IIIZw35 ( [Baan, Haschick, fc Henkel 1989| ). The 
positions of OH masers in several Galactic star-forming regions, e.g., W3(0H) ( [BloemhoT^ 



Reid, fc Moran 1992[ ), W75N ( [Baart et al. TOBQ and G5.89-0.39 QWood 1993[ ), have 



been observed to coincide with shock fronts in molecular outflows. The abundance of 
OH can be significantly enhanced in a post-shock medium (see, e.g., [Elitzur fc De Jong 
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1978; Hollenbach fc McKee 1989| ). However, unlike most extragalactic OH masers, those 



in Galactic star-forming regions typically have 1667 to 1665-MHz line ratios less than 
one ( [Uaswell, Haynes, fc Goss 198U| ). Mirabel, Rodriguez, & Ruiz (1989) have observed 
extended (pc-scale), weak OH maser emission in a number of star-forming regions that 
exhibits line ratios closer to those seen in extragalactic masers, and suggest that the two 
phenomena may be related; if so, then a shock origin for the observed maser structure would 



appear less likely. Detailed modelling ([Randell et al. 19951) also suggests that the brightest 



extragalactic OH masers may arise in pc-scale clouds with H2 densities on the order of 
lO^cm"^. We note that images of OH absorption toward the center of our Galaxy ( [Boyce 



Gohen 1994| ) reveal a number of "wideline clouds" that exhibit velocity widths in excess 



of lOOkms"^, and intrinsic velocity gradients. However, the typical size of these clouds is 
rather larger (~ 100 pc), and their velocity gradients at least an order of magnitude smaller, 
than in the maser clumps in IIIZw35. 

Although the continuous structure inferred by Montgomery & Cohen (1992) was 
probably a result of spectral blending, it is still possible that the two maser clumps lie in a 
rotating disk. We cannot rule out the presence of emission in the region between the two 
clumps, given that 50% of the total maser flux density was resolved. The velocity separation 
of the maser peaks could be indicative of rotation about a central mass of 2 x 10'' M©. A 
difficulty with this model is that the line profiles of the two clumps overlap in velocity by 
more than 50kms^^, which would imply superimposed non-rotational velocities comparable 
to the disk rotational velocity. Furthermore, the velocity gradients we observe are more 
than an order of magnitude greater than what we would expect from differential rotation 
alone. More elaborate models, such as a warped disk, or emission on the surface of a thick 
torus, cannot be ruled out by this observation, but are difficult to quantify. It is possible, 
despite the similar angular and velocity structure of the clumps, that they are physically 
unconnected structures, e.g., counter-rotating, gravitationally bound molecular clouds, or 
shocks associated with a molecular outflow. 



5. Conclusion 

We have used the VLBA to image the OH maser in the nucleus of IIIZw35. We detect 
maser emission in two clumps, with a projected separation of 50 pc along a position angle 
of —10°. The northern clump is extended, and exhibits a west-to-east projected velocity 
gradient of ~ llkms^^pc"^. The brightest maser emission arises in the southern clump, 
which is mostly compact (< 3pc); however, faint emission from the line wings suggests a 
velocity gradient of similar magnitude to that in the northern clump, but in the opposite 
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direction. We detected no 18 cm continuum emission in the nucleus, which places a lower 
limit on the apparent maser optical depth, ItohI > 3.4, assuming the maser is unsaturated 
and is amplifying continuum radiation. This is the highest apparent gain yet determined 
for an extragalactic OH maser. The observed flux ratios of the 1667 to 1665-MHz OH 
maser lines are consistent with optical depths of this magnitude or greater. Variations of 
the line ratio with angular position and velocity suggest that physical conditions are not 
uniform within the maser clumps. We did not find evidence for a 100-pc rotating, edge-on 
molecular structure, as reported by Montgomery & Cohen (1992); however, since 50% of 
the total maser flux was resolved, there is likely an extended or clumpy component of 
maser emission in IIIZw35 that does not appear in our images. Velocity gradients within 
the clumps may indicate rotation on smaller scales (3-10 pc). The observed gradients could 
be produced by two counter-rotating, gravitationally bound molecular clouds, each with 
an average density of roughly 7000 MqPC^^. Alternatively, the maser clumps could be 
segments of molecular shocks viewed edge-on. 
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Fig. 1. — Distribution of 1667-MHz OH maser emission in the nucleus of IIIZw35. The 
angular position offsets are referenced to the maser peak. Features greater than 5a (lOmJy) 
in the spectral images are indicated by circles; the diameter of each circle is proportional to 
the flux density of the feature. The heliocentric velocity of the spectral peak in each clump is 
indicated (optical definition; f (Heliocentric, Optical) = f(LSR, Radio) + 235kms~^ at fsys)- 
The directions and ranges of the velocity gradients observed in each clump are indicated with 
arrows. The synthesized beam is 25 x 12 mas with a position angle of —27°. Inset A plot 
of velocity vs. R.A. offset in the northern maser clump. The solid line corresponds to the 
best-fit velocity gradient (over the range 8290 < v < 8340 km s~^) of 5.6±0.4kms~"^mas~^, 
or 11 ± 1 kms~^ pc~^. 

Fig. 2. — Expanded views of the OH maser clumps: (a) The northern clump; (b) The 
southern clump. The diameter of each circle is proportional to the flux density of the feature. 
Selected features are labelled with their heliocentric velocities. The flUed squares indicate 
the positions of the 1665-MHz emission peaks; their hehocentric velocities are labelled in 
itahcs. 

Fig. 3. — The 1667-MHz OH maser spectrum of IIIZw35. The heavy solid curve is a single- 
dish spectrum taken from Staveley-Smith et al. (1987). The two solid curves show the total 
imaged power as a function of heliocentric velocity in each maser clump; the dotted line 
is the sum of both clumps. Arrows mark the velocities of the 1665-MHz emission peaks. 
The marker at 8289 km s~^ indicates the systemic velocity. Note that maser emission in 
single-dish spectra extends to velocities beyond those plotted here. 
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